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We derive an upper bound on the ground state energy of the three-dimensional (3D) repulsive 
Hubbard model on the cubic lattice agreeing in the low density limit with the known asymptotic 
expression of the ground state energy of the dilute Fermi gas in the continuum. As a corollary, we 
prove an old conjecture on the low density behavior of the 3D Hubbard model, i.e., that the total 
spin of the ground state vanishes as the density goes to zero. 

I. INTRODUCTION 

Recent developments in the theory of low density Bose and Fermi gases made it possible to 
verify old conjectures on the leading asymptotics for the ground state energy of dilute gases of 
bosons or fermions in the continuum, interacting with positive short range potentials. While 
the heuristing argument suggesting that the ground state energy of the 3D hard-core Bose gas 
is proportional to the scattering length a of the potential goes back to Lenz [Le] , the first ideas 
in the direction of proving that Lenz's formula is correct in the low density limit are due to 
Dyson [D], who first established an asymptotically correct upper bound and a rigorous (but 14 
times too small) lower bound for the hard core Bose gas in 3 dimensions. An asymptotically 
correct lower bound was proven much more recently by Lieb and Yngvason [LY]. Their work 
inspired much of the recent developments in the rigorous theory of low density quantum many 
body systems, see [LSSY] for a comprehensive review of the subject till 2005. In particular, a 
result that we would like to mention, strictly related to the problem studied in this paper, is 
the proof in [LSS] that the ground state energy per unit volume of the 3D Fermi gas in the 
continuum with short range repulsive interaction (and scattering length a > 0) is given, in the 
low density limit pa^ 0, by: 

e(PT,Pi) = i^l(^^')'^'ipT+pT) + I^S^PTPi (1-1) 

where p^^i axe the densities of spin up and spin down particles and m is their mass. Moreover 
p = Pi + Pl and o{ap^) is a suitable function of the total density p and of the scattering length 
a vanishing faster than ap^ in the limit pa^ ~+ 0. 

It is very natural to ask whether a formula similar to (1.1) is valid for a dilute Fermi gas with 
short range repulsive interaction on the lattice and, in particular, for the most popular model for 
correlated electrons in condensed matter physics: the Hubbard model. The Hubbard model is 
the simplest possible lattice model of interacting electrons displaying many "real world" features 
and in the last 40 years it has been subject of intense research efforts. Nonetheless, even its 
qualitative behavior in 2 or 3 dimensions is far from clear and there are very few rigorous results 
available in the literature. A survey of known results and open problems in the Hubbard model 
can be found in [Li] and in [T] . 

In the present paper we shall derive an upper bound for the ground state energy of the 3D 
repulsive Hubbard model on the cubic lattice with the same asymptotic behavior as (1.1). As 
a corollary we shall prove one of the open problems posed by Elliott Lieb in his review article 
on the Hubbard model (see [Li, Problem 3]). More precisely, we shall prove the following old 
conjecture on the low density behavior of the 3D Hubbard model. 

Proposition. Let Smax ~ ^totl"^ be the maximum spin a system of Ntot electrons can 
achieve and let S be the spin of the ground state of the 3D repulsive Hubbard model in a cubic 
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hox K clT? in presence of N tot electrons (or the maximum such spin in case of degeneracy). 
Then 

lim lim S/Smax = Q, (1.2) 

P^O |A|->oo 

where the thermodynamic limit is taken keeping the total density p = Ntot/\M fi-xed. 

Remark. A sketch of the proof of this claim aheady appeared in [BLT]: in this note we 
provide all the details of the necessary computations as well as an explicit bound on the rate 
of convergence (see Corollary 1 below). Our strategy imitates the one in [LSS]. 

The paper is organized as follows. In the next two subsections we shall introduce the model, 
introduce the notion of scattering length and state the main results, i.e., the upper bound on 
the ground state energy and an explicit bound on the rate of convergence of S/Smax to in 
(1.2). In Sec. II and in the two Appendices we shall give the proof. 



A. The model 



Given a cubic lattice A of lattice spacing ro, the Hamiltonian of the Hubbard model on A for 
N spin-up particles and M spin-down particles can be written as: 

H = -Ax -Ay + Uvxy (1.3) 

where: 

1) X — {xi, . . . jXff) and Y = (j/i, . . . ,yM) are the coordinates of the spin-up and spin-down 
particles, respectively; 

2) Ax = Eti Ax„ Ay = A,^. and AJ{x) = r^' E..|.'-.|=.o(/(^') - /(^)); 

3) vxY = Eili Ejli ^xi,vj; 

4) f/ > 0; 

5) H acts on the space of functions antisymmetric in the X and in the Y coordinates separately 
and vanishing outside the box A (Dirichlet boundary conditions). 

Remark. In this note we restrict for simplicity to the case of a nearest neighbor hopping 
and a delta interaction, however the analysis below can be generalized to cases with differ- 
ent hopping terms and different short range interactions (not necessarily zero - or finite - range). 

We want to obtain an upper bound for the ground state energy that is asymptotically correct, 
at the lowest order, as pa"^ — > 0, where p = {N + M)/\A\ and a is the scattering length of the 
potential. The latter can be conveniently defined in terms of the solution to the zero energy 
scattering equation 

-A^ip{x) + ^Sa^^oVix) = (1.4) 
subject to the boundary condition lim|a,|_>oo (fi{x) = 1. The solution is 

if{x) = l-4TT- . (1.5) 

ro4,|<,,-i (27rro-i)3 2E?=i(l-cosfciro) 
where the coefficient a has the interpretation of scattering length and is given by 

«^« = ^°Z7;^' ^=2y,,,.(2^2EL(l-cosfc.)- 

Note that lim|^|^(^(l — (^(a;))|a;| = a (this means that at large distances ^p{x) looks very much 
like the scattering solution in the continuum, i.e. if^ix) ~ 1 — a/\x\ at large distances) and that 
Stto <UrQ (this is the analogue of the inequality of Spruch and Rosenberg [SR] in the lattice 
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case). Another important remark is that, given any simply connected domain containing the 
origin, the "flux" of the discrete derivative of (f{x) across the boundary of f2 is independent of 
fl and equal to 47ra: 

(dQ) 

{<p{x') - ^{x)) = Ana (1.7) 

<x,x'> 

where "^^^^ x'y sum over the bonds connecting nearest neighbor sites with x ^ Q and 

x' gQ'^. This simply follows by the remark that A(f{x) =Q,^x ^ 0, and by discrete "integration 
by parts" . 

B. Main results 

We are now ready to state our main result. 

Theorem 1. Fix = N /\h\, pi = M l\K\ and p = p^ + pi, and let Eo{N, M, A) denote the 
ground state energy of H with the appropriate antisymmetry in each of the N, M coordinate 
variables. Then, for small pa^, 

Tn^o{N,M,A) <eo{p^,Pi) + 8TTap^pi+ap^s{pa^), (1.8) 

|A|^cx) |A| 

where eQ{pi,pi) is the specific ground state energy of the free Fermi gas on the lattice, i.e., of 
(1.3) with U = 0, andO < e{pa^) < const {p'^^^af/^. 

Remarks. 

1) The theorem is valid for any repulsion strength U >0, including the limiting case U = +oo 
of infinite repulsion. 

2) If prg ^ 1, the specific ground state energy of the free Fermi gas on the lattice can be written 

as 

eo(pT,Pi) = I ipf + pT) + const ry/' 

Then, as long as a/ro ^ \/ p^/^a, in the r.h.s. of (1.8) we can replace eo(pt)Pi) by 
I (67r^)^^'^ {p^^^ + p\^^) S'lid still have an error term that is much smaller than p^a. In this case 
the upper bound (1.8) looks precisely the same as (1.1). 

3) It would be nice to establish that the r.h.s. of (1.8) is the correct low density behavior of 
the ground state energy of the 3D Hubbard model. In order to prove this we should provide a 
lower bound with the same asymptotic behavior as the r.h.s. of (1.8). The natural idea would 
be to proceed as in the continuum case [LSS], that is by exploiting Dyson's idea of replacing 
the "hard" interaction potential by a "soft" one, at the expense of using up some kinetic 
energy. Of course, in order to get the correct 0-th order contribution in the lower bound, we 
need to use at least part of the kinetic energy to "fill the Fermi sea": so technically one of 
the main steps in the proof of the lower bomid in [LSS] is the proof of a "Dyson Lemma" in 
presence of an infrared cutoff, allowing for a replacement of the hard interaction by a soft one, 
at the expense only of the high momentum part of the kinetic energy. We would expect that 
this result is actually independent of the presence or absence of an underlying lattice struc- 
ture: however the proof of the "Dyson Lemma" with infrared cutoff in [LSS] uses in a crucial 
way rotational invariance of the problem and it is an open problem to adapt it to the lattice case. 

The result of the Theorem above, combined with the remark that the first term in the 
r.h.s. of (1.8) provides an obvious lower bound to the ground state energy, implies that 
\Eo{N,M,k) - e'^"^\n,M,K)\ < const|A|ap2 and this in turns implies that, if we fix the 
total density p and minimize the energy over the possible choices of pf^j^, we find that at low 
density the absolute ground state satisfies jpt — Px| < const p\/ p^/^a This implies that the 
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total spin S of the ground state satisfies the following. 

Corollary 1. Let S be the total spin in the absolute ground state of model (1.3) and let 
Smax = {N + M)/2. Then in the low density limit 

lim S/Sjnax < const -\/ p^/^a (1.9) 

|A|^oo 

where the thermodynamic limit is taken keeping the total density p = {N + M)/\A\ fixed. 

Remark. It is natural to ask whether there exists some number p,. > such that 
lim|A|^oo '5'/<S'moa: = for all p < Pc (see Problem 4 in [Li]). Note that Corollary 1 does 
not exclude this possibility. Note also that proving or disproving this possibility requires nec- 
essarily some non pertiirbative argument: any approximate computation of the ground state 
energy can only improve the error term in the r.h.s. of (1.9) but will never establish the exact 

value of S/Smax- 



II. THE UPPER BOUND 

In this section we shall assume a/ro > (5~^(p^/^a)^/^, with 5 a constant to be chosen be- 
low. In this case it is enough to prove the upper bound (1.8) with eo(/9t,/9j,) replaced by 

I (67r^) (/9| + P| ), sec Remark (2) after the statement of the Theorem above. The weak 

coupling regime a/ro < 5~^{p^/^a)'^^'^ is much simpler and will be treated in Appendix B. 

It will be convenient to localize the particles into small boxes with Dirichlet boundary condi- 
tions. The number of particles in each box will be large for small p, but finite and independent 
of the size of the large container A. Let the side length of the small boxes be i. We then want 
to put n = p-\l!? spin-up particles into each box, and likewise m = piU? spin-down particles 
(here p^ = A''|A|~-^ and pi = M|A|~^). Since Pi,ii^ need not be an integer, however, we will 
choose 

n = p|£^+£| and m = pi£^ + ei, (2.1) 

with < £|,| < 1 chosen such that n and m are integers. We then really have too many 
particles, but this is legitimate for an upper bound, since the energy is certainly increasing with 
particle number. 

So, if Eo{N, M, A) is the ground state energy of (1.3) in the box O, we have 

lim ^EoiN,M,A) < ^Eo{n,m,Ai), (2.2) 

|A|^oo |A| 

where A^ is the cubic box of side i. Here we used that the interaction potential is zero range, so 
that different boxes of side £ are exactly decoupled. Note that actually the bound (2.2) is not 
only valid in the thermodynamic limit, but also for all finite cubic boxes A with side divisible 
hyi. 

We will now derive an upper bound on the ground state energy of n spin-up and m spin-down 
particles in a cubic box of side length i, for general n, m and £. We take as a trial state the 
function 

*(X, Y) = Dn{X)Dm{Y)Gn{X)Gm{Y)F{X, Y), (2.3) 

where Z)„(X) denotes the Slater determinant of the first n eigenfunctions of the Laplacian in 
a cubic box of side length £, with Dirichlet boundary conditions; note that, if 0a (a;) are the 
eigenfunctions of the single-particle Laplacian in a cubic box of side length £, we choose their 
normalization in such a way that J2x''~o^a{^)^0{^) — ^a,0. Moreover, 



l<i<j<n 



(2.4) 
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with < g{x) < 1, having the property that g(x) = for < s and g{x) = 1 for \x\ > 2,s, for 
some s to be chosen later. We can assume that for any pair of nearest neighbor points x and 
x' we have \g{x') — g{x)\ < const ros~^ for some constant independent of s. Finally, 

n m 

F{X,Y) = l[l[f{xi-yj), (2.5) 

i=i j=i 

where, given a simply connected domain fl C 'Z'^ containing the origin, f{x) = 1 if x ^ fi. 
Inside we choose f[x) to be the solution to the zero-energy scattering equation with boundary 
conditions f{x) — ip(x) / {(p) i^^^ the boundary 9il of the domain (here ip{x) is given by (1.5), 
dQ = {x e f2 : dist(a:, fi'^) = 1} and {(p)gQ = \dil\^^ 'l2xedn fi^))- We shall make the following 
explicit choice for the domain: O = Bn fl Z^, where Bn = {x S : ip{x) < 1 — a/R}. Note 
that, if i? ^ ro, Bj^ is approximately a ball of radius R. Moreover {(p)qq = l~a/R+0{arQ/ R^) 
and, for any x G dft, (p{x) = {'f)gQ + O{aro/ R^). We assume 6~^ro < R < s/5, with 6 the 
same constant as in the condition a/ro > 5~^{p^/^a)'^/^ (to be chosen below). 
By the variational principle, 

£o(n,m,A,)<M^), (2.6) 

with 

= (*| - Axl*) + (*| - Arl*) + U^vxy]^)- 

(here, for any operator i, is defined as = Ejc.y ^o^"'^'"''*(^> ^) ^ ^) 

- note the presence of the factor Tq^"'*''"-'). We first evaluate {^\ — A^l^')- By definition it is 
equal to 

^ E E rf''^'"^D^{YfGm{YfDr,{X)Gn{X)F{X, Y) ■ (2.7) 

'"O i=l X,Y 

■ Yl [Dn{X)Gn{X)F{X,Y)-Dn{Xl)Gn{Xl)F{Xl,Y) 

x'^:\x'^—Xi\=rQ 

where, if X — {xi, . . . ,Xi, . . . ,x„j}, X^ is given by X[ = {xi, . . . ,x[,. . . ,x„j}. The r.h.s. of 
this equation can be written as 

^ rf''+"'^D^{YfGm{YfGn{XfF{X, YfDn{X){-Ax)Dn{X) 

X,Y 

+iEE^o^"^'"^ E D^{YfG^{YfD^{X)D^{X'^. (2.8) 

i=l X,Y xr:\x'.-Xi\=ro 

■ Gn{X)F{X, Y) [Gn{X)F{X, Y) - Gn{Xl)F{Xl, Y) 

The first line is simply E^{n,£){^\^), where E^{n,i) is the sum of the lowest n eigenval- 
ues of the Dirichlct Laplacian in the box of side C (note that these eigenvalues arc ciqual to 
2r(J'^ ^j_;^(l — cos/cjro), with ki positive integer multiples of n/i). An explicit computation 
shows that 

i?°(n,£) < ^(67r2)2/3!^ (l + const + const n2/3(ro/^)2) (2.9) 
The second line, if we symmetrize over X, X-, can be rewritten as 



1 " 

;iE E E r^^"+-)D„(y)2G„(r)2D„(x)D„(xo • (2.10) 

Gn{X)F{X, Y) - Gn{Xl)F{Xl, Y)] ' 



i=l X\xi,Y <Xi,x'^> 
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where X]<rE x'> ^^'^ o^'^'" ^^"^ nearest neighbor bonds in A^. By Cauchy-Schwarz, we can 
bound the last expression from above by 



1 " 

° i=l <Xi,xJ> 

■ 'Gn{X)F{X, Y) - Gn{X'^)F{X[, Y)] ' 
We now use the Schwarz inequality to deduce (for some e > to be chosen later) 

'Gn{X)F{X, Y) - Gn{X'^)F{X[, Y)] ' 

< (1 + e)G„(X02 [f{X, Y) - F{Xi Y)] ' 
+(1 + e-^)F{X, Y f \Gn{X) - G„(XOl ' 



(2.11) 



(2.12) 



Proceeding in the same way for the kinetic energy of the F-particles, we thus get the upper 
bound 



with: 



(*|lf|*)< [£;°(n,^)+£;°(m,£)](*|*) + (l+£)/2+(l + e-')/3, (2.13) 



h = Y.'-f^^^^Dn{XfD^{Yf 



X,Y 



=1 x'^:\x'^—Xi\=rQ 



F{Xi,Y) - FiX,Y) 



+G^{Xr\Yl E GrniY!)' 

J=i yy-\y'j-y3\='^o 
+UGniXfGmiYfvxYF{X,Yf I 



F{X,Y') - F{X,Y) 



(2.14) 



and 



^3 = J2^o"^^^D^iYfDr,{XrFiX,Yf 

X,Y 



Gm{Yf\VxG{X)\^ + Gn{Xf\VYG{Y)\' 



(2.15) 



where |VxG(X)|2 = Er=i |V^G(X)|2 and |ViG(X)|2 = i E^=± |VrG(X)|2. Moreover, for 
any fimction g{x), the tj-gradicnt of g is defined as V'^.g(.T) = \g{x + '^t'oc^) — 5(2;)] 

with ei the coordinate versor in the ^-th direction. A similar definition is valid for |VFG(y)p. 
The positivity of Uvxy has been used here. Note that E^{n,£) and E^{m,£) can be bounded 
as in (2.9). We shall now bound I2 and is, when divided by (^'|^'), separately. 

Let us first derive an upper bound on I2. We are going to need the following lemma [LSS]. 

Lemma 1. Let Dn{X) denote a Slater determinant of n linearly independent func- 
tions (t>a{x)- For a given function h{x) of one variable, let $(-X') he the function ^{X) = 
Dn{X) nr=i f^i^i)' '^'^^ ^ denote the nx n matrix 



M^^ = J2rlr.{x)Mx)\h{x)\' 



(2.16) 
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Then the norm of <I> is given by {^\^) = detM. Moreover, for 1 < k < n, the k-particle 

densities of $ are given by 

^)(iT5) S rr-"MX}P= (2.17) 
1 

= fc! n (a;i A • • • A a;fe jM"^ • • • (g) M"^ jxi A • • • A a;fe) , 



w/iere |a;) denotes the n-dimensional vector with components 4>a{x), 1 < a < and \xi A 

■■■Axk) stands for the Slater determinant X](t(~1)'^|2^o-(i)) ® ' ' ' ® ^(^(fc))' denoting 

permutations. Finally, if ^^{X) = Dn{X)k{xi) Ilj^j f^i^j) foi~ some function k{x), then 

n 

^($P^> = (detM) (Tr[™-i]), (2.18) 

where Tr[-] denotes the trace, and K is the nx n matrix 

Ka0 = Y.''o^-{^)M^)\H^W ■ (2-19) 

X 

Using Gni^i) ^ 1) we infer from this lemma that, for any fixed F, 



FiXi,Y) - F{X,Y)i2 



+ ^GniXfvxYF{X,Yy^ 

< E^o"^«(^)'||VxF(X,y)|2 + ^t,.;,y^^(X,y)2| (2.20) 
= Tr{KYMy'}Y,rl^F>r,{XfF{X,Yf 

X 

The n X n matrices -fCr and My are given by (2.16) and (2.19), with (l)a{x) being the lowest n 
Dirichlet eigenfunctions of —A, and with h{x) = JJ- f{x — yj) and 

\k{x)\' = |V/i(x)p + ^Y,Sx,y,l[fix - y,)\ 

3 3 

respectively (here \Vh{x)\'^ = {'^rQ)~^Ylix'-\x'-x\=ro ~ h{x')\'^, see definition after (2.15)). 

Since Ky is a positive definite matrix, we have the bound TrKyMy^ < \\My^\\TtKy, where 
II • II denotes the (spectral) matrix norm. To calculate TrKy, and to bound ||My^||, we can 
assume that all the y^'s are separated by at least a distance s, because the summand in (2.14) 
vanishes otherwise. 

Since s > 5R by assumption, we have in this case |A;(a;)|^ = ~ Vj) with 

ax -y) = |V/(a; -y)\^ + ^6x,yf{x - yf . (2.21) 
Hence, if p^i^) denotes the one-particle density of we have 

m m 

TrKy = EE^oP^(x)$(:r - yj) Y.Pn * ^(%)' (2-22) 
where * denotes convolution. In order to bound ||My^^||, we use the following: 
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Lemma 2. Assume that \yi — yj] > s > 5R for all i ^ j. Then if tq/R, n ^ and n{ro/t)^ 
are sufficiently small 

111 - My II < const (^+n^/^-^\ . (2.23) 



Proof. Let q{x) = 1 — Hj ~ Vj)'^ ^ 0- Then, for any n-dimensional vector |6) with 
components ba, 

{b\l-MY\b)=Y,rlq{x)\j2ba<f>a{x) 

X a 

Hence, the question about the largest eigenvalue of 1 — My translates into the question of how 
large the average potential energy for the potential q{x) can be for functions such as J2a ^a^Paix) 
whose kinetic energy is bounded above by (const.) n^/^^~^, i.e., the Fermi energy for n particles 
(under the assumption that n » 1 and utq <S £'^). 

Let Qj denote the cube of side s/2 centered at yj. Note that all these cubes are non- 
overlapping by assumption. Also, since s > 5i?, q{x) = if .x is outside all the cubes (recall 
that by definition - see the lines following (2.5) - f{x) is identically 1 outside a region Br 
of radius R[l + 0[{ro/R)'^)], for some k > 0). For a given function (j){x), let denote the 
average of (f>{x) in the cube Qj. Moreover, let ri{x) = (j){x) — By the Cauchy-Schwarz 
inequality (a + 6)^ < 2(a^ + 6^), we get the bound 

^ qixMx)]' < 2 ^ q{xMx)f + 2|</.,f Yl l{x). (2.24) 

Note that |(/>jp < 8s~'^ ^^^g^. ro|?ii(x)p, again by the Cauchy-Schwarz inequality. Moreover, 
since s > R, 

J2 rlq{x) = ^ rlil - f{xf) < const aR'. 

To obtain the last inequality, we used that if a; £ Br then f{x) = 'p{x)/{tp)g^ > ip{x), with 
(fi{x) defined in (1.5). 

Note that ri{x) is a function whose average over the cube Qj is zero. In other words, it is 
orthogonal to the constant function in Qj. Hence, using the fact that q{x) < 1: 

5^ r'oq{xMx)\' < J2 r'o\vix)\' (2.25) 

xeQj xeQj 



2(1 — cos(27rros )) 



where we used that 2r(7^(l — cos(27rro.s~'^)) is the second eigenvalue of the discrete Neumann 
Laplacian in the cube of side s/2 and mesh tq. In this last expression we can replace r]{x) by 
4>{x), of course, since they only differ by a constant. Summing over all the cubes Qj (and using 
that q{x) = outside the cubes), we thus obtain that, for any function <j){x), 



^r^g(x)|0(x)|2 < const 



In the case in question, the kinetic energy of <p{x) is bounded by const n^/^£~^. This finishes 
the proof of the lemma. ■ 
Since < My < 1 as a matrix, this lemma implies that 

l-||l-My|| -^"^ 1 -const [ai?Vs3 + „2/3 (5/^)2]' (2-26) 
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provided the denominator is positive. By inserting (2.22) and (2.26) into (2.20), we see that, 
for fixed Y with \yi — yj \ > s for all i ^ j, 

Y,rl^Dn{Xf^VxF{X,Y)\'' + ^vxYF{X,Yf^ 

n 

<AnY.Pn* ^{Vj) E rl-Dr.{XfF{X, Yf . (2.27) 

j=l X 

To be able later to compare this expression (2.27) with (\E'|4'), we want to put Gn{XY back 
into the integrand. For this purpose we need the following lemma, which compares the integrals 
with and without the factor G„(X)^. 

Lemma 3. For any fixed Y , if and n{ro/tf are sufficiently small 
J2 rl^D„{XfF{X,YfGn{Xf 

X 

>J2ro''D„{XfF{X,Yf (l - const n^/^WMy^fis/efy (2.28) 

X 

Proof Since g{x) = 1 for \x\ > 2s, we have 

n 

Gn{Xf>l-J2e{2s-\xi-Xj\). (2.29) 

Here 9 denotes the Heaviside step function, i.e., 9{t) = 1 for f > and 9{t) = for t < 0. To 
evaluate the sum involving the second term in (2.29), we need the two-particle density of the 

state Dn{X)F{X,Y) for each fixed Y. By Lemma 1 above, and the fact that f{x) < 1, this 

density, when appropriately normalized, is boimded from above by ||My "'^ ppj?'*'^'' (x, .t'), where 

Pn'''^\x,x') denotes the two-particle density of the determinantal state Dn{X). In particular, 
by explicit computation one finds that, if n 1 and nr^ <C this latter density satisfies the 
bound 

p°'(2)(a;,a;') < const \x - x'\''{n/ff/^ (2.30) 

for some constant independent of n and £. Hence we arrive at (2.28). ■ 
Let 



Bn=(l- const n^/^Al{s/£f^ ' , 



assuming that the term in parenthesis is positive. Applying Lemma 3 to inequality (2.27), we 
arrive at 



Y,rl^-+^^Gr^{XfD,,{Xf 

XY 

n 

3{n+m) D 



\yxF{X,Y)\^ + ^vxYF{X,Yf 



Dm{Y)Gm{Y) 



< AnB^ J2 E 4 Pn * ^{yj)Dm{YfDn{XfF{X, YfGm{YfGn{Xf . (2.31) 

J = l XY 

Now we cannot bound * ^(y) independently of y by simply using the supremum of p^{x), 
since this number will be strictly bigger than n/i'^ , even in the thermodynamic limit. Instead, 
we repeat the above argument for the Y integration. We use |Gm(5^)| < 1, the F-analogues of 
Lemma 1 and then Lemma 3 to put Gm(F)^ back in. Here, it is important to note that now 
the separated by at least a distance s > 5R. In this way we obtain 



X,Y 



\\/xF{X,Y)\^ + ^vxYF{X,Yf 



DmiY)G„,{Y) 



< Ar,Br,B„, J2 rf''^"'^Dm{YfDr,{XfF{X, YfGm{YfGr.{Xf TrKxM^' (2.32) 



XY 
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The matrix Mx is the same as before, with Y replaced by X (and n replaced by m, of course), 

and Kx is the m x m matrix 

(KxU = ^rlMvYMy) n - "^i^^Pn * ^iv) ■ 

V i 

Using |/(a:;)| < 1 and ||M^^|| < Am, which fohows from Lemma 2 and the fact that the 
separated at least by a distance s, we get the bound 



TrKxM^^ < AmTrKx < A„ ^ r6pD(a;)pD (y)e(x - y) . 

A computation (see Appendix A) shows that 

X! 'oC(2;) < 47ra(l + const a/R) . 



(2.33) 



(2.34) 



We then use this information to bound the last sum in (2.33), by using Schwarz's inequality: 

1/2 / V 1/2 



< f E ^IPni^m^ - y)] (e ^lp^n.{yn{x - y)\ 

\x,v / \x,y ) 



1/2 



1/2 



E^op!?(^)M E-op°(2/)M E-o?(-) 



(2.35) 



1/2 



^ \ 1/2 

< \^rlpl{xf]^ (^E-oP°(2/f j 47ra(l + const |). 
For the square of ^{x)-, by an explicit computation we find 

2 

^rlp^ixf < ^ (l+ const n-i/3 + const n2/3(ro/£)^) . (2.36) 

X 

The same holds with n replaced by m. Eq. (2.32) thus implies the upper bound 



E^o^"+'"^G„(X)2l?„(X)2 



XY 



\VxF{X,Y)\' + ^vxYF{X,Yf 



Dm{Y)Gm{Y) 



, , , Ananm . . „ „ 



(2.37) 



1 + const + + + („ + m)2/3(ro/£)2) 



The same bound holds, of course, with X and Y interchanged. We therefore have the upper 
bound 



(2.38) 



1 + const + « + m i/3 + (n + m)2/3(ro/^)2) 



It remains to bound the term 73. Using \g{x)\ < 1 we have that 



|VxG„(X)|2 < 1 E 



E E \'^''9{xi-xj)\' 



+ E E E \^^9{xi-Xj)-V'^9ixi-Xk)\ 



(2.39) 



11 



where the w-gradient was defined after (2.15). Now, by Lemma 1, the appropriately normal- 
ized fc-partiele densities of Dn{X)F{X,Y) are bounded above by ||M-^^||'^pJ?'*-'^\ where pn''''^^ 
denotes tlie fc-particle density of In particular, pn'^^^ satisfies the bound (2.30), and 

Pn'^^^ satisfies 

p^''^^\x,x',x") < const {n/i'^f 

for some constant independent of n and i. Remember that, by the definition of g (see lines 
following (2.4)), if x and x' are two neighbor points, g{x') — g{x) is zero for |a:| > 2s and 
otherwise \g{x') — g{x)\ < const roS~^. Using this, we obtain from (2.39), for any fixed Y, 

Y,rl''D4XfF{X,Yf\WxGn{X)\'' (2.40) 

X 

< const (\\My'fn'/^s/£f + \\My'\\Ms/if) Y,rl"Dr,{XfF{X,Yf . 

X 

Finally, to get a bound on I^, we proceed as above, using (2.26) (and the fact that the j/j's are 
separated by a distance s) and Lemma 3 to put G„(X)^ back into the integral. Note, however, 
that it is enough to bound A„ and Bn by constants in this term. Assuming that n{s/i)^ is 
small, the second term in the parenthesis in (2.40) is negligible compared to the first term. The 
same bound applies to the case where X and Y are interchanged, and hence we obtain 

I3 const (n^/^ + ^. (2.41) 

Collecting all the error terms obtained in Eqs. (2.9), (2.38) and (2.41) and inserting them 
into (2.6) and (2.13), we obtain 



Eo{n,m,Ae) < 



(1 + Cn-V3 + (7m-V3 + C{n + mf/^{ro/lf^ 



£3 



+87ra— 5- l + e + C 



^ + (n + m)2/3(s/£)2 



1 



+ 



R ni/3 rnV 
Cs (n + mf - 



L_ + (n + mf'\s/lf 



(2.42) 



£3 



(n + m)2/3(s/£)2 



for some constant C > 0. In Ineq. (2.42) we have assumed smallncss of all the error terms, 
i.e., that the terms in square brackets are small. This condition will be fulfilled, at low density, 
with our choice of R, s, n, m and i below. 

The optimal choice of e in (2.42) is given by = const {n + m)^^^s^/{i'^anm}. Inserting this 
value for e we infer from (2.42) 



Eo{n,m,Ai) < 



(1 + Cn-V3 + (7m-V3 + C{n + mf'\r<y/lf) 



+87ra 



nm 

IF 



aR^ 

c3 



+ (n + m)2/3(s/£)2 



(2.43) 



+C(n + m)^/3 



,3/2^1/2 



Eq. (2.44) is our final bound on the energy EQ{n,m,£). To apply this result in (2.2) we have 
to insert the values (2.1) for n and to. Recall that \n — p^i^\ < 1 and |to — Pi£^\ < 1. We are 
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then still free to choose R, s and i. We choose 

R = a(apV3)-V9 ^ ^ ^ ^ ^ ^ p-V3(apV3)-iV9_ 

Note that with these choices the condition a/ro > S~^{p^/^ay^^ imphes ro/R < 6. We choose 
S to be so small that Lemma 2 is valid. Inserting these values into (2.44) we thus obtain, for 
small p, 

^Eoin, m, Ke) < ^(eTr^)^/^ + ^5/3^ ^ Sivap^p^ + const ap^ {ap^l^f'\ 

In combination with Eq. (2.2), this finishes the proof of the upper bound in the case a/ro > 
5~^{p^/^a)'^^^ . The opposite case (that is much simpler) is treated in Appendix B. 
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APPENDIX A: PROOF OF (2.34) 



In this Appendix we want to show that if £^{x) is defined by (2.21) then 'I'qS^{x) < 47ra(l + 
const a/-R). Note that "integrating by parts" we can rewrite the summation as: 



r^/(x)(-A,)/(x) 



(A.1) 



Note also that, by definition, if x G O, then f{x) coincides with <p(a;)/((^)gQ, while if a; ^ f2 
then /(x) = 1. Then, by the definition of (p{x), wc see that in the summation in (A.1) all 
terms with x "well inside" Q (i.e. with x such that dist(a;, fi"^) > 2ro) cancel out with ■^/(O)^, 
and all terms with x "well outside" fl (i.e. with x such that dist(a;, O) > 2'ro)) are identically 
zero. So we are left with a boundary term, that is a summation over the x's at a distance 
ro from Q or from 0°. Let us recall that dfl = {x € fl : dist(a;, fi'^) = ro} and let us define 
dQ" = {x gCI'^ : dist(a;,f2) = ro}. The r.h.s. of (A.1) can be rewritten as 



^ r3/(x)(-A)/(x) + J2 rlf{x'){-A)f{x') 

(^) r / /\ 1 (^') 



(A.2) 



E 

xedn 



ro 



E 



an 



E 



if{x) 



an 



where Xli^eon= is the sum over the points x' G 90'^ at a distance ro from x (and similarly for 



X^iean)- "^^^ second line in (A.2) can still be rewritten as 



E 

xean 



ro 



an 



ix) 

E 

x'ean= 



an 



(V')an 



E ro^ix') - ^{x)] (A.3) 



<x,x'> 



where 'YT^xx'y '^^ ^^"^ sum over the nearest neighbor pairs < x,x' > with x 6 dVl and x' G 
dff^. Recall that, if i? ^ ro, then for any x € dil and any x' € 90° we have (p{x),(p{x') = 
i'P)dn + 0{aro/R^)- Then the first term in (A.3) can be bounded above by a constant times 



(i?^/rQ)ro(aro/-R^)^ = a{aro/R'^) < a{a/R). Moreover, note that the second term in (A.3) 
is proportional to the (discrete) flux of the discrete derivative of across the "surface" of 
and the latter is equal to 47ra, see (1.7). As a conclusion, the second term in (A.3) is equal to 
47ra/((/3)gQ. Using that (v')an = 1 ~ a/R + 0{aro/R^), see lines following Eq. (2.5), (2.34) is 
proven. 
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APPENDIX B: THE WEAK COUPLING REGIME 

In this Appendix we prove the main Theorem in the case that o/ro < 5~'^{p^^'^a)'^/^ . In this 
case wc do not locahzc particles and we simply choose as trial function the ground state of the 
free Fermi gas: '^{X,Y) = Di^{X)Dm(X), where Djq{X) denotes the Slater determinant of 
the first N eigenfunctions of the Laplacian in the cubic box A with (say) periodic boundary 
conditions (a similar definition is valid for DMiX))- We assume that Di^{X) and Dm(^) 

arc normalized in such a way that (*|*) = Ejc.y ^o^^'^^^I-Dw(^)P|£'m(>")P = 1- By the 
variational principle i;o(iV,M, A) < = 4^=°^ W A) + E.li Ej=i« |*). 

Since the specific energy corresponding to the term Eq{N, M, A)|(7=o is by definition eo(pi-, pj,), 
we are left with bounding 

N M 
Y i=l X j=l 

An apphcation of Lemma 1 shows that this expression is equal to J2x'''qPn{^)pm{x), where 
pN is the 1-particle density of Dn{X) and Pm is the 1-particlc density of Dn{Y). In our case 
Pn{x) = pt and pm{x) = py. Then we get C/E^i Ej-i" |*) = \h.\Urlp^pi. % (1-6) 
we have that Ur^ = 87ra(l + ^Uro). As a conclusion: 

lim i^i;o(7V,M,A) <eo(pT,Pi) + 87rapTP4(l+7[/r^) (B.2) 

lAHoo |A| 

Since Utq = 87ra/ro(l +consta/ro), we have that ^Ur^ < const (p^/^o)^/^ and the proof is 
concluded. 
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